function res=fx(t,R)
q=1000;
pM=101.325*10^3;
pv=2.318*10^3;
sigma=0.0728;
u=1002*10^(-6);
ye=1.4;
A=0.4;
R0=0.26*10^(-3);
res=zeros(2,1); %res保存函数表达式，t对应时间，R对应半径和速度，R1、R2—>R,R'
res(1)=R(2);    %高阶化为低阶,R1=R,R1'=R2
res(2)=(-3/2*R(2)^2-1/q*(pM-pv-(A*pM-pv+2*sigma/R0)*(R0/R(1))^(3*ye))-(4*u/q)*(R(2)/R(1))-2*sigma/(q*R(1)))/R(1); %R2'=R1''=R''，关键表达式
end